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Abstract 

Our recent study, based on a 3D numerical simulation, suggested that the stochastic excitation 
of solar p-modes has a non gaussian nature and a non gaussian model was proposed which 
results in a good agreement between the computed rate P at which p-modes are excited and 
the solar seismic observational constraints. 

In the present work, we study some consequences of this non-gaussian model on the rate 
P at which p-modes are excited in intermediate mass stars (1M© < M < 2Mq). The non 
gaussian model changes substantially the spectrum of P at high frequency compared with 
the commonly assumed gaussian model. The largest effects on P(y) are found for stars in 
the mass range ~ 1.5 — 1.6 M®. They are found large enough that observatins in solar-like 
oscillating stars observed with space based seismology missions (e.g. COROT) will be able to 
confirm or reject the gaussian nature of stochastic excitation. 

1 Introduction 

Excitation of solar-type oscillations results from turbulent movements in the upper convective 
zone of intermediate mass stars (1M Q < M < 2M Q ). The rate P{v) at which a given mode 
with frequency v is excited crucially depends on the time averaged and dynamic properties of the 
turbulent medium (Samadi 2001). 

Different models for stochastic excitation have been proposed by several authors (e.g. Goldreich 
k Keeley 1977; Balmforth 1992; Goldreich et al. 1994; Samadi & Goupil 2001). The model 
of Samadi & Goupil (2001, Paper I hereafter, see also Samadi (2001)) takes into account the 
frequency dependence Xk{v) of the correlation product of the turbulent velocity field. This last 
quantity characterizes the dynamic properties of the turbulence. This excitation model offers 
then the advantage to test different hypothesis concerning the dynamic properties of the turbulent 
medium. 

A gaussian function is usally assumed for Xk(y) ( e -g- Goldreich & Keeley 1977). However 
recently, using 3D numerical simulation of the upper part of the solar convective zone, Samadi et al. 
(2002) have shown that the gaussian model indeed does not correctly model Xk{v) in the frequency 
range where the acoustic energy injected into the solar p-modes is important (v ~ 2 — 4 mHz). 
One must consider an additional non-gaussian component for Xkiy) to reproduce its behavior. 
Computed values of P obtained with this non-gaussian component reproduce better the solar 
seismic observations. 
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In the present work we study the consequences of this non-gaussian component for the stochas- 
tic excitation of p-modes in intermediate mass stars. Sect. 1 describes the stochastic excitation 
model we consider and sumarizes Samadi et al. (2002) 's results which establish the non gaussian 
character of the stochastic excitation. In Sect. 2 we compute P for different stellar models with 
masses in the mass range 1M Q < M < 2M . In the computation we assume two hypotheses for 
Xk(y) : the gaussian function and a non gaussian function as inferred in Samadi et al. (2002) from 
a 3D simulation. 



2 The non gaussian character of the stochastic excitation 

The theoretical expression in Paper I for the excitation rate P can be written in a schematic form 
as : 

P(v) oc J dm J dxdr £ . < SS> (r,r).£ (2.1) 

where r and r are related to the local turbulence : r is the distance between two given points in 
the turbulent medium and r is the time correlation between the two distant points (r is associated 
with the eddy lifetime which is typically much shorter than the oscillations lifetime) . In Eq. ( {0] ) , 
£ r denotes the radial mode eigenfunction and < S S > the correlation product of the excitation 
source versus r and r. The excitation source (S) originates from the turbulent Reynolds stress 
and the advection of the turbulent entropy fluctuations by the turbulent motions. < S S > is 
expressed in terms of the turbulent kinetic energy spectrum E(k, v) and the spectrum of the 
entropy fluctuations E s {k,v) where k is the wavenumber of a given turbulent element. Further 
E(k, v) is split into a spatial component E{k) and a frequency component Xkiy) according to Stein 
(1967) as E(k, v) = E(k) Xk{ v ) The same decomposition is assumed for E s {k,v). 

In the calculation of P{v), a gaussian shape is usually assumed for Xk{v) (e.g. Goldreich & 
Keeley 1977). This assumption is equivalent to suppose that two distant points in the turbulent 
medium are uncorrelated. Such Gaussian function (GF hereafter) has the form 

xTiv) = ^e~W-> 2 , (2.2) 
Vk V 77 

where is the line-width at half maximum. 

With the help of a 3D numerical simulation of the upper part of the solar convective zone, it is 
found in Samadi et al. (2002) that the GF does not reproduced correctly the frequency dependence 
of Xk inferred from the 3D simulation. In contrat, Xk is better modeled by the Gaussian plus an 
Exponential function (GEF hereafter) as: 

- \ (-^e-^ + ^-e-^A . (2.3) 

2 V^feV 71 " 2ly k 1 



3 Consequences for solar-like oscillating stars 

We compute several stellar models with masses M = 1., 1.2, 1.5, 1.65, 1.85, 2. M© and with the 
same age (650 My). Calculations are carried out with the CESAM evolutionary code in which 
the convective heat flux is computed according to the classical mixing-length theory by Bohm- 
Vitense (1958). The eigenfunctions are obtained from the adiabatic FILOU pulsation code by 
Tran Minh & Leon (1995). The input physics are similar to those of Samadi et al. (2001) except 
for the atmosphere which is computed with the Eddington classical gray atmosphere for sake of 
simplicity. 



We compute P for each the stellar model according to Eq. (2.1). In the calculation, we assume 



successively for Xk(v) : the gaussian model (GF, Eq. 2.2) and the non gaussian model (GEF, 



Eq. 2.3). A scaling factor is introduced such that for the solar model the maximum in P matches 
the solar seismic constraints by Chaplin et al. (1998). Results are presented in Fig. [l]. 

For all stellar models, the largest differences in P between calculations performed with the 
GF and those performed with the GEF are obtained for v > f max where i/ max is the frequency at 
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Fig. 1. The excitation rate P at which stellar p-modes are excited is computed for several stellar models 
with masses M = 1.2, 1.5, 1.65, 1.85 Mq. Two different analytical functions are assumed for Xfc(^) : the 
GF (dashed curve, Eq. 2.2) and the GEF (solid curve, Eq. 2.3). Vertical error bars AP are estimated in 
the manner of Samadi et al. (2001). The assumed accuracy in the measurement of the mode damping 
rates (r/) and the mode relative luminosity fluctuations (SL/L) are based on the expected performances 
of the COROT mission (Baglin et al. 1998) (i.e. Ar]/2n = 0.1 fj,Uz and ASL/L = 0.6 ppm). Values of r] 
are calculated in the manner of Houdek et al. (1999). 



which P is maximum. This feature is a direct a consequence of the difference in the frequency 
dependence of the GF and the GEF. Indeed, the GEF decreases with v slower than the GF does 
for v > V}.. Moreover ^ max scales approximatively as where is the characteristic frequency 
associated with the largest energy-bearing eddy located at the top of superadiabatic (this is the 
region where the excitation is the largest). The eddies who excite the most the modes have their 
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characteristic frequency ~ v^. Consequently for v > f m ax a larger amount of acoustic energy 
is injected with the GEF than the GF. 

The largest effects in P of the dynamic model adopted for \k are obtained for the M = 1.5 Mq 
and M = 1.65 Mq stellar models : Let AP be the difference in P between calculations performed 
with the GF and those performed with the GEF. AP is found larger for these two models than for 
the other stellar models. Moreover AP takes large values over a larger range in v/v c where v c is 
the cut-off frequency (y c delimits the extend of frequency domain where modes have a significant 
amplitude). These features are due to the fact that the characteristic turbulent Mach number, 
M t , is the highest for stars with mass ~ 1.5 — 1.6 Mq. Indeed, AP max , the maximum in AP, 
scales approximatively as v\ A\k where Axk = ||Xfe EF ~ Xfe F ||i/=p fc - As Ax^ scales as D^ 1 and Dk 
scales approximatively as M t v c , AP max increases then as M 2 . Moreover, let Av be the frequency 
domain where || AP— AP max | < e AP max where e is arbitrarily fixed to a value <C 1. In other words 
Av defines a domain where AP is "large". One shows that d 2 Axk/dv 2 scales as D^ 3 and then 
d 2 APjdv 2 oc ^,7 2 AP max . As a consequence ||AP - AP max || w | D^ 2 Av 2 AP max . The domain 
Avjv c then scales as y/euk/vc oc M t . 



4 Conclusion 



The present work illustrates that a non gaussian model of stochastic excitation substantially 
changes the shape and the magnitude of the calculated excitation spectrum P of the oscillations in 
intermediate mass star more massive than the Sun. Differences in P between computation carried 
out with the non gaussian excitation model and this carried out with the common gaussian model 
are found larger for the frequency modes above the frequency ^ max at which P peaks. The largest 
effects are expected for stars in the mass range ~ 1.5 — 1.6 Mq. They are found large enough to 
be discriminate from space based experiences with performances equivalent or higher than those 
of the COROT intrument (Baglin et al. 1998). 
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